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Abstract —This paper presents an omnidirectional spatial and 
temporal 3-dimensionaI statistical channel model for 28 GHz 
dense urban non-line of sight environments. The channel model 
is developed from 28 GHz ultrawideband propagation measure¬ 
ments obtained with a 400 megachips per second broadband slid¬ 
ing correlator channel sounder and highly directional, steerable 
horn antennas in New York City. A 3GPP-like statistical channel 
model that is easy to implement in software or hardware is 
developed from measured power delay profiles and a synthesized 
method for providing absolnte propagation delays recovered from 
3-D ray-tracing, as well as measnred angle of departure and angle 
of arrival power spectra. The extracted statistics are used to 
implement a MATLAB-based statistical simulator that generates 
3-D millimeter-wave temporal and spatial channel coefficients 
that reproduce realistic impulse responses of measured urban 
channels. The methods and model presented here can be used 
for millimeter-wave system-wide simulations, and air interface 
design and capacity analyses. 

Index Terms —28 GHz millimeter-wave propagation; channel 
modeling; multipath; time cluster; spatial lobe; 3-D ray-tracing. 

I. Introduction 

Millimeter-waves (mmWave) are a viable solution for allevi¬ 
ating tbe spectrum shortage below 6 GHz, thereby motivating 
many recent mmWave outdoor propagation measurements 
designed to understand the distance-dependent propagation 
path loss, and temporal and spatial channel characteristics of 
many different types of environments m-ia. The mmWave 
spectrum contains a massive amount of raw bandwidth and 
will deliver multi-gigabit per second data rates for backhaul 
and fronthaul applications, and to mobile handsets in order to 
meet the expected 10,000x demand in broadband data over the 
next 10 years ID 1^ . 

MmWave statistical spatial channel models (SSCMs) do 
not yet exist, but are required to estimate channel param¬ 
eters such as temporal multipath delays, multipath powers, 
and multipath angle of arrival (AOA) and angle of depar¬ 
ture (AOD) information. Both directional and omnidirectional 
channel models are needed based on real-world measurements. 
Further, SSCMs are needed to carry out link-level and system- 
level simulations for analyzing system performance required 
for designing next-generation radio-systems. New channel 
modeling frameworks, modulation schemes and corresponding 
key requirements are currently being considered to address 
5G network planning |2l-||9l, and future directional on-chip 
antennas m require that 3-D statistical channel models be 


used at mmWave bands. This paper presents the world’s first 
3-D SSCM based on wideband New York City measurements 
at 28 GHz. 

Previous results obtained from the extensive New York City 
propagation database yielded directional and omnidirectional 
path loss models in dense urban line of sight (LOS) and non¬ 
line of sight (NLOS) environments El, important temporal 
and spatial channel parameters, and distance-dependent path 
loss models at 28 GHz and 73 GHz based on measurements 
and ray-tracing results 121 na. Initial MIMO network simu¬ 
lations were carried out in 1T3]| using a 2-dimensional (2-D) 
wideband mmWave statistical simulator developed from 28 
GHz wideband propagation measurements 12], and showed or¬ 
ders of magnitude increase in data rates as compared to current 
3G and 4G LTE using spatial multiplexing and beamforming 
gains at the base station for both LOS and NLOS dense urban 
environments. 

Statistical channel modeling methods have thus far focused 
on extracting models from measured power azimuth spectra 
and by modeling the elevation dimension using 3-dimensional 
(3-D) ray-tracing predictions in order to make up for the lack 
of measured elevation data El or to estimate 3-D angles in 
the absence of directional measurements due to the use of 
quasi-omnidirectional antennas na. In this paper, we present 
a 3-D statistical spatial and temporal channel model for the 
urban NLOS environment based on 28 GHz ultrawideband 
propagation measurements that used 3-D antenna positioning 
and directional antennas in New York City, where AOA 
elevation characteristics at the receiver have been extracted 
from measured power angular spectra (obtained with field 
measurements without the use of ray-tracing). This work 
extends our 2-D SSCM presented in 121 to a true 3-D model 
that includes the AOA of arriving multipath, as well as AOD, 
and mutipath time statistics. 

H. 28 GHz Propagation Measurements 

In 2012, an ultrawideband propagation measurement cam¬ 
paign was carried out at 75 transmitter (TX) - receiver (RX) 
locations in New York City to investigate 28 GHz wideband 
channels using a 400 megachips per second (Mcps) broadband 
sliding correlator channel sounder and 24.5 dBi (10.9° half¬ 
power beamwidth (HPBW)) highly directional steerable horn 
antennas at the TX and RX, in dense urban LOS and NLOS 
environments m. Over 4,000 power delay profiles (PDPs) 
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Fig. 1: Superimposed PDF of two individual received PDFs, where 
each PDP comes from a different AOA at the same RX location. 
The multipath signals from Angle 1 arrived before those of Angle 
2 (i.e. multipath arriving at different times from two distinct lobes). 
The absolute propagation times were found using ray-tracing, thus 
allowing alignment with absolute timing of multipath signals origi¬ 
nally measured at the RX, independent of AOAs (see ©)• 



Fig. 2: A 3-dimensional view of the downtown Manhattan area using 
the MATLAB-based 3-D ray-tracer. The rays which leave the TX 
and successfully arrive at the RX are shown in red, and represent 
multipath signal paths. The TX was located on the rooftop of the 
Coles Sports Center 7 m above ground (yellow star), and the RX 
was located 113 m away, 1.5 m above ground (black circle) (!]. 


were collected at unique azimuth and elevation pointing an¬ 
gles at both the TX and RX to properly describe mm Wave 
wideband multipath channels over an 800 MHz RF null-to- 
null bandwidth in the time and spatial domains. Additional 
details pertaining to the 28 GHz measurement campaign and 
hardware equipment used can be found in IT] ||3 ifTTI ifTSl . 

HI. Synthesizing 3-D Absolute Timing 
Omnidirectional Power Delay Proeiles 

The collected directional PDPs did not make use of absolute 
timing synchronization between the TX and RX over the 
various sweeps, and therefore provided received power over 
excess, and not absolute, time delay. As illustrated in Fig. [T^, 
two typical measured excess time delay PDPs are shown as 
measured at two distinct AOA azimuth angles. The RX used 
the strongest arriving multipath component to trigger and 
establish the relative f = 0 ns time marker for all recorded 
PDPs, as illustrated in Case 1 of Fig. [TJ) where both PDPs 
from different angles are shown to start at f = 0 ns, and 
thus could not capture absolute arrival time of power. Case 2 
illustrates the desired result, showing the two PDPs measured 
along an absolute propagation time axis, where f = 0 ns 
corresponds to the energy leaving the TX antenna. In Case 2, 
the recorded excess time delay PDPs have been shifted in time 
appropriately, accounting for the propagation distance travelled 
of the first arriving multipath in each PDP, enabling accurate 
characterization of received power in time at the RX over all 
measured angles. 

We developed a MATLAB-based 3-D ray-tracer to predict 
the propagation time delays of the first arriving multipath 
components at the strongest measured AOAs for each TX- 
RX measured location. Ray-tracing techniques have previ¬ 
ously been shown to provide accurate time and amplitude 
predictions l4] IfT^ IfTSlI . Fig. shows a typical ray-traced 
measured location where viable propagation paths between 
TX and RX are shown in red. The ray-tracing results showed 
strong spatial correlation with the measurements when com¬ 


paring the strongest measured and predicted directions of 
departure and arrival (within ± two antenna beamwidths). The 
predicted propagation distances were paired with the closest 
strongest measured AOAs, and we superimposed each PDP 
recorded at the strongest AOAs in azimuth and elevation on 
an absolute propagation time axis by appropriately shifting and 
summing each PDP in time using the ray-tracing absolute time 
predictions. This method was performed over all measured 
RX locations, and allowed us to synthesize omnidirectional 
absolute timing PDPs as would have been measured with a 
quasi-omnidirectional isotropic antenna, yielding one unique 
omnidirectional 3-D power delay profile at each TX-RX loca¬ 
tion combination. Our measurement approach enabled high 
gain antennas to measure the channel over large distances 
in a directional mode, and the orthogonal beam patterns in 
adjacent beam pointing angles allow the orthogonal PDPs to 
be summed in time and space to form a 3-D omnidirectional 
model. 

IV. 3-D Statistical Channel Parameters 

Our statistical channel model recreates omnidirectional 
channel impulse responses h^,~^ where t denotes abso¬ 
lute propagation time delay, © = {9, (j))TX represents the vec¬ 
tor of AOD azimuth and elevation angles, and ^ = {S, (j))RX 
is the vector of AOA azimuth and elevation angles, as defined 
in Eqs. and (|2l. Previous work modeled the channel 

impulse response as a function of time uni, azimuth angle 
of arrival m HD, and AOD/AOA elevation and azimuth 
dimensions ll20l . Note that ll20l includes AOD and AOA 
elevation and azimuth information for in-building, indoor- 
to-outdoor, and outdoor-to-indoor scenarios, but not for the 
outdoor urban microcellular environment as considered in this 
paper. Here, we generalize the channel impulse response as 
a function of time, as well as a function of AOD and AOA 
azimuth/elevation angles, allowing realistic simulations of di¬ 
rectional transmissions at both the TX and RX in azimuth and 

































elevation dimensions. The channel impulse response between 
TX and RX is written as: 

Lapp Lao a N M„ 
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where L^od and L^qa are the number of AOD and AOA 
spatial lobes (defined in m ED, and in Section [3), re¬ 
spectively; N and M„ are the number of time clusters and 
the number of intra-cluster subpaths in the n* time cluster 
(as defined in Section lIV-Bl i. respectively; am,n,h,i 2 is the 
amplitude of the m* subpath component belonging to the 
n* time cluster, departing from AOD lobe li, and arriving 
at AOA lobe I 2 ', ‘Pm,n, and tm,n are the phase and the 
propagation time of arrival of the to* subpatl^ component 
belonging to the n* time cluster, respectively; 0;^ and 
are the azimuth/elevation AODs and AOAs of lobes li and I 2 , 
respectively. In our channel model, each multipath component 
is assigned one joint AOD-AOA lobe combination, per the 
time cluster definition given in Section IIV-BI 

Our statistical channel model also produces the joint AOD- 
AOA power spectra in the azimuth and elevation domains in 
3-D, based on our 28 GHz New York City field measurements 
that only used a single TX pointing elevation angle at a 10° 
downtilt and three RX elevation planes of 0°, ±20° (note: 
measurements at 73 GHz provided multiple AOD and AOA 
elevation angles, and may be extrapolated and used here). 
The spatial distribution of power is obtained by taking the 
magnitude squared of and integrating over the 

time dimension, as shown in (|2]i and Q: 

p{4,4)= r ( 2 ) 
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Note that P(^,4) in Eq. Q are the total received powers 
(obtained by integrating the PDF over time) assigned to the 
lobe AODs and AOAs. 

A. 28 GHz Omnidirectional NLOS Path Loss Model 

The 28 GHz omnidirectional NLOS path loss model was 
recovered by summing the received powers measured at each 
and every azimuth and elevation unique pointing angle com¬ 
bination to recover the total received omnidirectional power 
at each TX-RX location ifTTl . This procedure is valid since 
adjacent angular beamwidths are orthogonal to each other, 
and phases of arriving multipath components with different 
angular propagation paths can be assumed identically and 
independently distributed (i.i.d.) and uniform between 0 and 
such that powers in adjacent beam angles can be added. 
After removing antenna gains and carefully removing double 
counts occuring from the TX and RX azimuth sweeps, we 
recovered the corresponding path loss at all measured NLOS 


locations, and extracted the path loss exponent and shadow 
factor using the 1 m close-in free space reference distance 
path loss model ||4l ifTTl : 

PLNLOs{d)[dB] = 61.4± Sdlogipjd) +Xa, d > 1 m (4) 

where Xa- is the lognormal random variable with 0 dB mean 
and shadow factor a = 9.7 dB, and 61.4 dB is 28 GHz free 
space path loss at 1 m. 

B. Cluster and Lobe Statistics 

The temporal and spatial components of our SSCM are 
modeled by a time cluster and spatial lobe, respectively, and 
faithfully reproduce omnidirectional PDFs and power azimuth 
spectra in. Time clusters model a group of multipath com¬ 
ponents travelling closely in time over all directions, and can 
represent one or more spatial or angular directions within the 
same time epoch. Spatial lobes represent a small contiguous 
span of angles at the RX (TX) where energy arrives (departs) 
over a small azimuth and elevation angular spread. In our 
SSCM, multiple time clusters can arrive in one spatial lobe, 
and a time cluster can arrive over many spatial lobes within 
a small span of propagation time. Time clusters and spatial 
lobes statistics can be easily extracted from the propagation 
measurements to build a 3GPP-like statistical channel model, 
including simple extensions to the current 3GPP and WINNER 
models that account for cluster subpath delays and power 
levels II 22 II . Time cluster and spatial lobe characteristics were 
illustrated in ||9l (see Figs. 3 and 4). The measured data 
suggests that multiple clusters in the time domain arrive up 
to several hundreds of nanoseconds in excess time delay for 
arbitrary pointing angles, observable due to our 24.5 dBi high 
gain antennas. We note that multipath components within 
a time cluster, i.e., intra-cluster subpath components, were 
successfully used to model the indoor office environment 
based on wideband measurements cii m. 

Key parameters that serve as inputs to the 3-D mmWave 
SSCM are referred to as primary statistics, and have been 
identified and defined in ||9l, and include the number of time 
clusters and the cluster power levels. Secondary statistics 
describe statistical outputs of the SSCM, and include the RMS 
delay spread and RMS angular spreads which reflect second- 
order statistics. Secondary statistics provide a means of testing 
the accuracy of a statistical channel model and simulator over 
a large ensemble of simulated outputs. 

C. Time Cluster Partitioning 

In this work, omnidirectional PDFs were partitioned in time 
based on a 25 ns minimum inter-cluster void interval, by 
assuming that multipath components fall within time clusters 
that are separated by at least 25 ns in time span. Walkways 
or streets between building facades are typically 8 m in 
width (roughly 25 ns in propagation delay). The 25 ns inter¬ 
cluster void interval allowed us to resolve measured multipath 
channel dynamics in a simple, yet powerful way, offering a 
scalable clustering algorithm that can be modified by changing 
the inter-cluster void interval for arbitrary time resolution. 






As the minimum inter-cluster void interval is increased, the 
number of time clusters in a PDP is expected to decrease, 
while the number of intra-cluster subpaths must consequently 
increase. In turn, the fewer time clusters must be allocated a 
larger portion of the total received power, while the greater 
number of intra-cluster subpaths will receive a lesser amount 
of the cluster power. While the minimum inter-cluster void 
interval heavily affects the outcome of the temporal channel 
parameter statistics, we note that the RMS delay spread is 
the only time parameter that remains unchanged for arbitrary 
inter-cluster void interval, making it a good but not the only 
proper indicator for comparing the simulated PDP outputs 
with the ensemble of measured PDPs. The 25 ns minimum 
inter-cluster void interval of our model is comparable to the 
best/maximum measured time resolution (20 ns) of a single 
multipath component in 3GPP and WINNER models. 

Fig. [3 shows the temporal cluster powers normalized to the 
total received power in the omnidirectional profiles, and the 
least-squares regression exponential model that reproduces the 
measured cluster powers. The cluster powers were obtained by 
partitioning omnidirectional PDPs based on a 25 ns minimum 
inter-cluster void interval, finding the area under each time 
cluster, and dividing by the total power (area under the 
PDP). The mean exponential curve is parameterized using two 
parameters, the average cluster power Pq in the first received 
cluster (i.e., the y-intercept at r = 0 ns), and the average 
cluster decay time T defined as the time required to reach 
37% (1/e) of Pq. It is worth noticing the measured large 
cluster power at r = 80 ns, containing close to 80% of the 
total received power, corresponding to large fluctuations in 
cluster powers. This phenomenon causes large delay spreads, 
and is typically modeled using a lognormal random variable, 
as discussed in Step 5 of Section |V] In Fig. 13 we estimated 
Pq = 0.883, and F = 49.4 ns. Similarly, Fig. |4] shows the 
intra-cluster subpath power levels (normalized to the total 
cluster powers), with Pq = 0.342 and 7 = 16.9 ns. The 
smaller subpath decay time physically means that intra-cluster 
subpaths decay faster than time clusters. 

D. 3-D Lobe Thresholding 

The 3-D spatial distribution of received power is used to 
extract 3-D directional spatial statistics, by defining a -10 dB 
lobe threshold below the maximum peak power in the 3- 
D power angular spectrum, where all contiguous segment 
powers above this threshold in both azimuth and elevation 
were considered to belong to one 3-D spatial lobe. Spatial 
thresholding was performed on 3-D AOA power spectra, and 

2- D AOD azimuth spectra, after linearly interpolating the 
directional measured powers (in linear units) to a 1 ° angular 
resolution in azimuth and elevation domains to enhance the 
angular resolution of our spatial statistics. 

V. Generating 3-D mmWave Channel Coefeicients 

A statistical channel model is now presented for generating 

3- D mmWave PDPs and spatial power spectra that accurately 
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Fig. 3: Temporal cluster powers normalized to omnidirectional 
total received power, over cluster excess delays, using a 25 ns 
minimum inter-cluster void interval. The superimposed least-squares 
regression exponential curve has an ave^e cluster decay constant 
of r = 49.4 ns, and a y-intercept of Po = 88.3% (see Step 7 
in Section IV-Al l. A large cluster power can be seen at r = 80 ns, 
which typically causes large delay spreads. 
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Fig. 4: Intra-cluster subpath powers normalized to time cluster 
powers, using a 25 ns minimum inter-cluster void interval. The 
superimposed least-squares regression exponential curve has an 
werage cluster decay constant of 7 = 16.9 ns, and a y-intercept of 
Po = 34.2% (see Step 8 in Section FV-Al l. 


reflect the statistics of the measurements over a large ensem¬ 
ble, valid for 28 GHz NFOS propagation with a noise floor 
of -100 dBm over an 800 MHz RF null-to-null bandwidth, 
and for a maximum system dynamic range of 178 dB m. 
The clustering and lobe thresholding methodologies used in 
our work effectively de-couple temporal from spatial statistics. 
Step 12 of our channel model bridges the temporal and spatial 
components of the SSCM by randomly assigning temporal 
subpath powers to spatial lobe AODs and AOAs, thereby re¬ 
coupling the time and space dimensions to provide an accurate 


















joint spatio-temporal SSCM. In the following steps, DU and 
DLN refer to the discrete uniform and discrete lognormal 
distributions, respectively, and the notation [cc] denotes the 
closest integer to x. Also, Steps 11 through 15 apply to both 
AOD and AOA spatial lobes. 

A. Step Procedure for Generating Channel Coefficients 

Step 1: Generate the T-R separation distance d ranging from 
60-200 m in NLOS (based on our field measurements, and 
could be modified with further measurements): 

d'^U niform{drmn = 60, d^nax = 200) (5) 

Note: To validate our simulation, we used the above distance 
ranges, but for standards work any distance less than 200 m 
is valid. 

Step 2: Generate the total received omnidirectional power 
(dBm) at the RX location using the 1 m close-in free space 
reference distance path loss model m^IB■ 

Pr[dBm]= Pt + Gt + Gr-PL[dB] (6) 

PLjvLOs[dB] = 61.4 + 34logio (d) + X> 7 , d > 1 m (7) 

where Pt is the transmit power in dBm, Gt and Gr are the 
TX and RX antenna gains in dBi, respectively, A = 0.0107 m, 
and ri = 3.4 is the path loss exponent for omnidirectional TX 
and RX antennas im. X(T is the lognormal random variable 
with 0 dB mean and shadow factor a = 9.7 dB. 

Step 3: Generate the number of time clusters N and the 
number of spatial AOD and AOA lobes (Laod, Lao a) ot the 
TX and RX locations, respectively: 

A^-D[/[l, 6 ] ( 8 ) 

Laod ~ min|Lmax,max|l,min{A, A^}|| (9) 

Laoa ~ min|Lmax,max|l,min{B, A^}|| (10) 

and: A ^ Poisson(iiAOD + 0.2) (11) 

B ~ Poisson(/ryioA + 0.1) (12) 

where haod = 1-6 and p,AOA = 1-7 are the mean number 
of AOD and AOA lobes observed in Manhattan, respectively, 
and Lmax = 5 is the maximum allowable number of lobes, for 
both AODs and AOAs. Work in ll2Tl found the mean number 
of AOA lobes to be 2.5 using a -20 dB threshold, whereas here 
we use a -10 dB threshold. The 28 GHz NLOS measurements 
found the maximum number of clusters Nmax = 5 based on 
measurements in ill], while at 73 GHz we found N^ax — 
6 , therefore we used N^ax — 6 for both frequency bands. 
Note that (Laod: Laoa) must always remain less than or 
equal to N, since the number of spatial lobes must be at most 
equal to the number of traveling time clusters in the channel. 
Also, the use of coin flipping was introduced in our previous 
work 191 to generate the pair (N, Laoa), to obtain a close fit 
between measured and statistical data. In this work, however, 
we use standard well-known distributions, without the use of 
coin flipping, to promote ease of use in simulated software. 


Step 4: Generate the number of cluster subpaths (SP) Mn in 
each time cluster: 


Mn^DU[l,30] , n = l,2,...N (13) 


At 28 GHz, the maximum and second to maximum number 
of cluster subpaths were 53 and 30, respectively, over all 
locations, while it was 30 at 73 GHz, so we choose to use 
30 for both frequency bands. 

Step 5: Generate the intra-cluster subpath excess delays pm,n' 


Pm,n(,Bi}}f) — \ ^ ipi 1 ) r 

L Dbb j 


l-hX 


(14) 


where Bf,b — 400 MHz is the baseband bandwidth of our 
transmitted PN sequence (and can be modihed for different 
baseband bandwidths), X is uniformly distributed between 0 
and 0.43, and m = 1, 2, ...Mn, n = 1,2, ...N. This step allows 
for a minimum subpath time interval of 2.5 ns, while reflecting 
our observations that the time intervals between intra-cluster 
subpaths tend to increase with time delay. The bounds on the 
uniform distribution for X will likely differ depending on the 
site-specific environment, and can be easily adjusted to fit field 
measurement observations. 

Step 6: Generate the cluster excess delays Tn (ns): 


< ^ Exp(/x^) (15) 

At„ = sort(T") - min(r") (16) 


Tn 


0 , 

Tn-1 + PM„_i,n-l + Ar„ -f 25, 


n = 1 

n = 2,..., N 


(17) 


where = 83 ns, and sort() in (fThl l orders the delay elements 
r" from smallest to largest. This step assures no temporal 
cluster overlap by using a 25 ns minimum inter-cluster void 
interval. 

Step 7: Generate the time cluster powers Pn (mW): 


P'^ = P^e-^ID^ (18) 

Pn = p, X Pn[mW] (19) 

Z^fe=l Tf. 

- A(0,3 dB),n = 1,2,...A (20) 

where T = 49.4 ns is the cluster decay time, Pg = 0.883 


is the average (normalized) cluster power in the first arriving 
time cluster, and is a lognormal random variable with 0 dB 
mean and cr = 3 dB. Eq. (fT9l l ensures that the sum of cluster 
powers adds up to the omnidirectional power P^, where Pg 
cancels out and can be used as a secondary statistic to validate 
the channel model. 

Step 8: Generate the cluster subpath powers Ilm,n (mW) : 


U'„,^n = ^oe-^10^ ( 21 ) 

H' 

Hm.n = X Pn[mW] (22) 

Z^k—l ^k,n 

Um,n N{0,6 dB) (23) 

m=l,2,...,M„ , n=l,2,...7V (24) 





where 7 = 16.9 ns is the subpath decay time, Hq = 0.342 
is the average subpath power (normalized to cluster powers) 
in the first intra-cluster subpath, and Um,n is a lognormal 
random variable with 0 dB mean and cr = 6 dB. For model 
validation, the minimum subpath power was set to -100 dBm. 
Eq. d22li ensures that the sum of subpath powers adds up to 
the cluster power. Note; our measurements have much greater 
temporal and spatial resolution than previous models. Intra¬ 
cluster power levels were observed to fall off exponentially 
over intra-cluster time delay, as shown in Fig. |4]and in ||9l . 
Step 9: Generate the cluster subpath phases iprn,n (rad) : 


V^m,n ~ 

- 4^1,n 4“ 27r f p^ n 

(25) 


(7(0, 27r) 

(26) 

m = 2, 

...,Mn,n = 1,2,...,N 

(27) 


where / = 28 x 10® Hz, and Pm,n are the intra-cluster subpath 
delays in s from Step 5, where / can be any carrier frequency. 
The subpath phases iprn,n are i.i.d and uniformly distributed 
between 0 and 2tt, as modeled in El El. 

Step 10: Recover absolute time delays tm,n (ns) of cluster 
subpaths using the T-R Separation distance: 


lm,n fo 4” '^n 4“ Pm^r, 


tn — — 


(28) 


where m = 1, 2, ...Mn, n = 1,2, ...TV, and c = 3 x 10® m/s 
is the speed of light in free space. 

Step 11 a: Generate the mean AOA and AOD azimuth angles 
0i(°) of the 3-D spatial lobes to avoid overlap of lobe angles: 


19U\0jnin, Omax] 

360(f - 1) 


i=l,2,...,L 
_ 360i 

^m.ax — Z 


(29) 

(30) 


Step 11 b: Generate the mean AOA and AOD elevation angles 
of the 3-D spatial lobes: 


f,r^[N{p,a)] , i = l,2,...,L. (31) 


where 5rs corresponds to the Kronecker delta. The cluster 
subpath [m, n), with power level Tlm,n, is assigned to lobe 
i only if Wm,n = *■ This step distributes subpath powers into 
the spatial domains based on measurements in ET\ . 

Step 13: Generate the AOA and AOD lobe azimuth and 
elevation spreads Ki (azimuth) and Hi (elevation): 

For AODs: K, - max(5°, [iV(30°, 16°)]), H, = 10° 

For AOAs: ~ DLN{32°, 18°), H, ~ max(5°, [N{3r, 11°)]) 

AOD elevation spreads are fixed at 10° based on our 28 GHz 
measurements that used a 10.9° antenna beamwidth. We have 
allowed for at most 10% lobe azimuth and elevation overlap 
in adjacent spatial lobes. 

Step 14: Generate the discretized lobe segment azimuth and 
elevation angles (Oij, fi^i) for lobe AODs and AOAs: 

~ 3- kj, j = 1,2, ...Ki, i = 1,2, ...L (34) 

4>i,i = 4>i + hi,l = 1, 2, ...Hi, *=1,2, ...L (35) 

where: {X, W) - DU[Q, 1],Y = 1-X,Z=1-W. 

\k, = -i^,...,-l,ll,l,...,i^, iV, odd 
\k, = -i^+Y,...,-l,0,l,..,k^-X, iV, even ^ ^ 


=-2^,...,-l,0,l,...,2f^, iJ, odd 

[hi = -I- W,..., -1,0,1,..., ^ - Z, Hi even 

This step discretizes the spatial lobes into 1° angular segments 
in both azimuth and elevation dimensions. 

Step 15: Generate the AOD and AOA lobe angular powers 
P{9ij,(j)i^i){mW) at each 1° angular segment: 

= R{xe,^j,X(t>^,i)P{e„ f,) (38) 

f + 1 ] 

R{A9ij,A(j)iy) = may:le V ^ ’ iQ J 

j = 1, 2,..., K,, 1 = 1,2,..., H,, i = l,2,..., L (40) 


Positive and negative values of fi indicate a direction above 
and below horizon, respectively. While our 28 GHz Manhattan 
measurements used a fixed 10° downtilt at the transmitter, and 
considered fixed AOA elevations planes of 0° and 4i20° at the 
receiver, mm Wave transmissions are expected to beamform in 
the strongest AOD and AOA elevation and azimuth directions 
as was emulated in our 73 GHz Manhattan measurements in 
Summer 2013 ||2|. Thus, we specify (/*, u) = (—4.9°, 4.5°) 
for AOD elevation angles, and {p,a) = (3.6°, 4.8°) for AOA 
elevation angles from our 73 GHz NLOS measurements. 

Step 12: Generate the AOD and AOA lobe powers P{9i, ff) by 
assigning subpath powers llm,n successively to the different 
AOD and AOA lobe angles (9i, fi): 

N M„ 

P(0„(/),) , i = l,2,...,L (32) 

n=l m—1 

Wm,n ~ DU[1, L\ and Srs = ^ ^ (33) 

l,r = s 


For AODs, ae, - TV(6.6°, 3.5°) and cr^, - TV(5°,0°). 

For AOAs, ag, - W(6°, 1°) and - A^(6°, 2°). 

B. Measured vs. Simulated Statistics using a MATLAB-Based 
Statistical Simulator 

A MATLAB-based statistical simulator using our 3-D 
SSCM generated a large ensemble (10,000) of mmWave PDFs, 
and AOD and AOA power spectra, where simulated and 
measured channel statistics were compared. Fig. |5] shows 
the cumulative distribution function (CDF) of the synthesized 
and simulated RMS delay spreads, showing relatively close 
agreement, with a median of 31 ns and 32 ns, respectively. 
The RMS delay spread CDF was skewed due to one large 
value of 222.4 ns, and so the median was chosen to reflect the 
empirical distribution trend. Fig. |6] shows the CDF of the RMS 
lobe azimuth and elevation spreads from the measured and 
simulated AOA power spectra, showing very good agreement 
to field measurements with identical measured and simulated 
means of 7° over both azimuth and elevation. 
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Fig. 5: 28 GHz NLOS synthesized and simulated omnidirectional 
RMS delay spreads. The median RMS delay spreads were 31 ns and 
32 ns, for the synthesized and simulated data sets, respectively. The 
close agreement between measured and simulated data validates our 
NLOS statistical channel model. 


28 GHz NLOS Lobe RMS Angular Spreads 



RMS Lobe Angular Spreads (°) 


Fig. 6: 28 GHz NLOS RMS lobe azimuth and elevation spreads, 
measured as compared to simulated, using a -10 dB lobe threshold. 
The simulated data is in good agreement to the measured RMS 
angular spreads, validating the spatial component of our NLOS 
statistical channel model. 


VI. Conclusion 

This paper presents the first comprehensive 3-D statistical 
spatial channel model for mmWave NLOS communication 
channels. The thousands of measured PDFs have allowed us 
to create a 3-D SSCM that recreates the measured channel 
statistics over a large ensemble of simulated channels, and 
can be extended to arbitrary bandwidths and antenna patterns 
for use in physical layer simulations, such as physical layer 
design, and 3-D beamforming and beamcombining simulations 
used in MIMO systems in NLOS environments. 
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